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Abstract 

The magnetic shielding a of ^He is studied. The complete relativistic corrections of order O(a^), leading 
QED corrections of order 0{a^ In a), and finite nuclear mass effects of order 0{'m/m^) are calculated with 
high numerical precision. The resulting theoretical predictions for a = 59.967 43(10) • 10~^ are the most 
accurate to date among all elements and support the use of '^He as a NMR standard. 
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I. INTRODUCTION 



The coupling of the nuclear magnetic moment /I to an external magnetic field B in closed shell 



atoms is shielded by electrons and thus is slightly smaller in comparison to the free nucleus Jl, 
This shielding is described by the dimensionless constant a 



m 



H = -fl- B{l-a). (1) 

For a specified atom cr is a function of the fine structure constant a and depends also on the 
electron-nucleus mass ratio. Since this ratio is very small, (for ^He it is about 1.8 ■ 10"^), a can 
be expanded in a power series in m/m^. We demonstrate in this work that the leading term in 
the mass ratio is numericaly significant, amounting for ^He to —3.7 ■ 10^"^ a, which is 42% of the 
relativistic correction. Regarding dependence on a, as long as the nuclear charge Z is sufficiently 
small, say Z < 10, the expansion in a is also well convergent. We have therefore for a a double 
series expansion 

a = a(a,—)=a' cr^^) + a' a^^^ + a^^) + a' —a^^''^ + ... (2) 

The first term of this expansion, a^'^\ is obtained from the Ramsey nonrelativistic theory llD of the 
magnetic shielding and for atomic systems takes the very simple form shown in Eq. (fTT)) . The 
derivation of the next coefficient a^"^^ was considered in a series of works by Vaara and collab- 
orators [|3|]. They expressed a^^^ in terms of the first, second, and third order expectation values 
of certain operators with the nonrelativistic wave function. Numerical evaluations of a^'^^ were 
performed for various elements, but somehow not for ^He, for example by Ruud et al. [4]. These 
calculations were not complete, in the sense that the authors omitted some terms which correspond 
to Q^, Qs, Qii from our Table HI all of them come from the Breit interaction (the second term in 
Eq. (|7])). These terms are small, but nevertheless important for the estimation of theoretical un- 
certainties. Moreover, the inclusion of the exact electron g-factor instead of the factor 2 by the 
authors of Ref. [jsl], in our opinion, is incorrect, as we explain in Sec. Ill, devoted to QED effects. 
In general, the relativistic correction cx^^-* can be obtained from the Breit-Pauli Hamiltonian in Eq. 
(|5]), in a very similar way to (t*^^^ In this work, apart from evaluation of the complete a^"^^ and 
(j(2,i)^ we present the calculation of the leading logarithmic QED correction cr^^\ which has a nu- 
merical value of about 10^^ cr. Finally we consider further improvement of theoretical predictions 
for ^He and the other light atomic systems. 
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II. RELATIVISTIC CORRECTION 



The relativistic correction a^^^ can be derived from the generalized Breit-Pauli Hamiltonian, 
which in addition to relativistic corrections includes coupling to the external magnetic field ||7|]. 
In our case, this field consists of the magnetic field Aj coming from the magnetic moment of the 
nucleus 

^1 r 

(3) 



r 1 ^ f 

^/ = ^ ^ X — 

Arc 



and of the homogenous external magnetic field Ae 

1 



Ae = - B X r . 
^ 2 



(4) 



This Breit-Pauli Hamiltonian in natural units with the external magnetic field A and with q being 

n 

the electron g-factor is DTD 



H 



BP 
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(7) 



where 71 = p — eA, and Hbp in the above includes dependence on the electron g-factor only for 
spin dependent terms. For the derivation of a^'^^ we set (7 = 2 and separate Hbp into parts: the 
leading interaction with the external field, no external field Hap] x a r., linear in the homoge- 
nous field SxHnplx linear in the nuclear magnetic field 67 Hnpl „, and bilinear in the 
homogenous and the nuclear magnetic fields S^^ ^^Hbp- 

o2 



Hbp = -Ae-Ai-^ J^iLa + 2sa)-B 
m 2m ^-^ 



+H^ 



bp\x^=Aj=q + ^AeHbp\ 



Ai=Q 



+ hHBp\ 



+ 5 



AeAjHbp 



(8) 
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Corrections to the energy, which are biUnear in magnetic fields, can be represented in terms of the 
leading contribution E^^^^ and the relativistic correction E^^\ 



m 



(9) 



7 7 1 
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m [E — Hy 



H 



BP 



+2 ( 5,Hbp 



1 



\Ae=Ai=0 

5xH 



Ai=o i^E - H)' ^'^I^B=o 

+ {hE,AHBp) . 



(10) 



Because of the spherical symmetry of closed shell atoms, perturbations due to — YlniLa + 
2sa) ■ B vanish, which is a significant simplification over molecular systems. Thus, in the leading 
order the shielding constant a takes the form 



while relativistic corrections are 

r(4) 



' = a 
2 
3 



W , ^(4) , ^(4) , (4) (4) 
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(16) 



1 

12 



r ■ n r ■ r2 



ri ■ r2 



(17) 



where r = ri — r2, (p* g*)*^^-' = /2+p^ q' /2 — 5'^ p ■ q/3, and 1/{E — H)' is the reduced Green 
function (the reference state is subtracted out). We have split these relativistic corrections into first 
order terms 0-3'*^ the second order terms with intermediate singlet ^S- a[^\ triplet ^S- a^J, singlet 
^P^- (72^, and triplet ^D- a^^c- These terms form a complete relativistic correction of order 0{a^) 
and their numerical calculations are described in Sec. V. 
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III. QED EFFECTS 



The next order correction O(a^) comes from QED effects. They contribute by Fi electromag- 
netic formfactors llisll . the magnetic susceptibility and the so called Bethe logarithms. The slope 
of -F{(0) and ^2(0) = [g — 2)/2 are known analytically at one-loop order [|l3n . However, F[{Q) 
is infrared divergent and this divergence cancels out with the uv divergence from the low-energy 
contribution in a similar way as for the Lamb shift in hydrogen [|l3ll . The contribution from ^2(0) 
is encoded by the g-factor in the Breit-Pauli Hamiltonian Hbp, Eqs. (|5l|7]). We note that the g- 
factor enters relativistic corrections with different coefficients, which is not in accordance with 

n 

Ref. [4]. The radiative corrections to the magnetic susceptibility have not yet been evaluated, 
but their calculation can be performed along the lines of Ref. [llSll . Finally the Bethe logarithmic 
contribution, as for the Lamb shift, is probably the most difficult part of the numerical evaluation 
and can be obtained probably only for simple systems. However, the total QED correction can 
easily be estimated on the basis of the leading logarithmic contribution, which is derived below. 

The leading logarithmic correction to a can be obtained in the same way as with the Lamb shift. 
One considers 2-electron self interaction in the magnetic field due to low-energy photons (m = 1) 



E, 



(2 7r)3 2fc 



,7r 



E-H-k 



(18) 



where (p is the eigenstate of the nonrelativistic Hamiltonian H in the magnetic field, with energy 
E. The logarithmic part is 



E, 



Log 



a 
a 

3^ 
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AnZa6{fa) + 2e^Bl + -{na,Va x 5,} 



(19) 



The resulting logarithmic contribution to the shielding constant in order is 



a 



(5) 



ln[(^«)-]((^ + ^)(^[^(n) + ^(r.)] 
ln[(Z«)-2] (5(fi) + 5(r2)> 



8Z 
20 

+ -— ma [ d{r 
9 



1 1 

(E-H)' Vn 



(20) 



where we have added a small second order In a term, which comes from the two-electron Lamb 
shift liO], and thus Eq. (|20|) is a complete logarithmic contribution. 
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IV. NUCLEAR MASS CORRECTIONS 



Effects coming from the finite nuclear mass are usually, if not always, neglected. Here we use 
the result derived in our previous work [|61], namely the leading correction a'^^'^) is given by 

^( > ^ 



(2,1) 



Ta [E - Hy I 3 Z gN 



(21) 



+ ^ X P2 + f2 X Pi) ^^3^ 



Era 



where Pn = - Y.aPa^ ^"^^ 



gN 



mN /i 1 



Zmp /in / 

where /^n is the nuclear magneton, / is the nuclear spin, and nip is the proton mass. The definition 
of the nuclear g-factor assumed here is different from the standard one by the use of the actual 
charge ejv and the mass mjy of the corresponding particle, namely the coupling of the spin to 
the magnetic field is — cat (7Ar/(2 rriAr) / ■ B. The numerical value of obtained from the known 
spin J = 1/2, nuclear charge Z = 2, magnetic moment /i = —2.127 625 2(1) fij^, and the mass 
ratio m^/rrip = 2.993 152 671 3(26) is presented in the caption of Table|Ill while the numerical 
evaluation of a^"^'^^ is performed in Sec. V. Our result is not in agreement with the work on Neronov 
and Barzakh in [Sfl, which based on earlier results of Hegstrom in [Qfl. 



V. NUMERICAL CALCULATIONS 



The numerical calculations are performed with the use of explicitly correlated exponential func- 
tions, which for the S-state have the form 



[ri,r2,r) 



Vile 



-air\-l3ir2-^ir 



± (ri ^ r2)], 



(23) 



where ctj, Pi and 7^ are generated randomly with conditions: 



Ai<ai< A2, A + 7i > £, 
Bi< /3i< B2, + 7i > e, 
Ci < 7i < C2, ai + /3i> e. 



(24) 



with e approximately equal to -y/2m (-E'hp+ — i^Hn)- In order to obtain a more accurate wave 
function, following Korobov [I11I1 . we use double set of the form (|24l) . Parameters Aj, Ci are 



determined by minimization of the nonrelativistic energy. The linear coefficients Vi in Eq. (1231) 
are obtained from a solution of the generalized eigenvalue problem with the length of the basis set 
Af = 100, 300, 600, 900, 1200, 1500 using extended precision arithmetic. As a result we obtain the 
following nonrelativistic energy in au 

Eo{l^So) = -2.903 724 377 034 119 593(5), (25) 



in agreement with the even more accurate result of Korobov nlu\ and of Drake in Ill2n . The 
calculation of matrix elements of the nonrelativistic Hamiltonian are performed with the use of a 
simple formula for the master integral: 

rfVi / £r2 = 7 — ^ r- r. (26) 



IGvr^ 7 J ri r2 r (a + + 7)(7 + a 

Integrals with any additional powers of in the numerator can be obtained by differentiation 
with respect to the corresponding parameter a, (3 or 7. Matrix elements of relativistic corrections 
involve inverse powers of ri,r2,r, and these can be obtained by integration with respect to the 
corresponding parameter. In fact, all matrix elements involved in relativistic, QED and the finite 
nuclear mass corrections can be expressed in terms of rational, logarithmic and dilogarithmic 
functions in a, (3, and 7. Considering numerical convergence, first order matrix elements can 
be calculated accurately using even short, i.e. M = 300 expansion of the nonrelativistic wave 
function. The calculation of the second order corrections Qe-is.is-ie is more complicated. The 
inversion of the operator E — H is performed in a basis set of even parity functions with / = 0, 1, 2 
of the form in Eq. (|23]) and 

— * 

0(ri,r2,r) 



J^^kTiX r2 [e-"'=-i-'3'=-2-7.'^ _ ^ r2)] , (27) 

k 

Vi, r2, r) = Y,^k [(ri 4 - r\ 5'V3) e-"^-'^^-^'^^^-^'^^ - (n ^ r^)] 

k 

■^^vi (r* - 5'^/?,) [e-^'^'i-^''^^-^''^ - (n ^ r2)] (28) 
I 

The values of parameters Ai, Bi and Ci are obtained by minimization of the appropriate functional, 
which is the symmetric second order matrix element with the operator standing on the right hand 
side of the corresponding Qi in Table I. This explicitly correlated exponential basis set allows us 
to obtain precise matrix elements of all Qi operators, and the numerical results are presented in 



Table I. Some of these matrix elements have already been presented in [|l2ll . and results in Table I 
are in agreement with them. 



TABLE I: Expectation values of operators entering a, all digits are significant, r = n — r2, and 1/{E — H) 
is the nonrelativistic Green function. 



^ + i 



1 „2 , J_ 2 
Pi ^2 ^2 



Qi 

Q2 

Qa 



(nxpi)^ , (f2Xp2)^ 



J(fi) +(5(f2) 



r-r2 
3 



Q6= [<5(n)+5(r2)](^(^ + ^) 

^9= <^(-"') (ft + ^) 

QlO = (pM X PI +plf2 X P2) ^ + ^) 



(r 2 X #1 + r 1 X p2 ) - r 1 X r 2 ^ • (pi + P2 ) J ( ^ + ^ ) 

n — (iiii (2) ]. /' 7 »i£i 7 '•2''2 I i:M!i±ii)^(2) 

Ql2 - ^) [^^-^^ + f3 ) 

Qi3 = [5in) - 5{r2)] ^ {3pI - 3p1 - 5 + 5 - ^) 

Qu = (pi +P2f 

Ql5= (Pl+P2f (^(^ + ^) 

gi6= (rixp2 + r2xpi)^(^ + ^) 



3.376 633 601 
33.677 743 

0.073 109 

3.620 859 
-3.435 251 
-3.025 857 
-118.140 232 
-0.208 449 
-0.112 964 
-0.041 690 

0.037 546 
-6.838 1(4) 
-39.921 269 

6.125 588 
-3.504 997 

0.078 743 



VI. RESULTS 



All corrections to the shielding constant a can be expressed in terms of Qi values and results 
are presented in Table IL All of them are accurate to all digits shown, nevertheless the uncertainty 
is different from 0, due to the neglect of the non logarithmic part of a^^'> which we estimate, on the 
basis of the helium Lamb shift, to be about 20%. The relativistic correction cr^^^ is relatively large, 
namely 10~^ of the nonrelativistic one, and is dominated by the second order contribution from the 
triplet 5" states. QED corrections are non negligible, 1% of the relativistic contribution, while the 
finite nuclear mass corrections are very significant, about 42% of the relativistic contribution and 
of the opposite sign. Except for neglected QED contributions, in our opinion, no other correction 
including the finite nuclear size, may alter the result at the 0. 1 ppb level. 
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TABLE II: Contribution to the shielding constant. Physical constants are taken from il4ll : a ^ = 
137.035 999 679(94), = -6.368 307 2,m^/m = 5 495.885 276 5(52). Uncertainty of fj is set to 

20% o f o-(^) contribution 

operator expectation value contribution to cr x 10^ 



a(^) = 


1 / ^ 
3 ^1 




1.125 544 534 


59 


936 


770 


(4) 


^ Qe — j2 Q7 - I Qs + ^ Qd 




3.340 6 








(A) 

41 = 


— 9- Q13 




27.870 3 










-g Qio - I Qii 




0.000 7 










-| Q12 




0.854 8 








(T3 — 


— 1 Q2 — 1 Qs — 3- <54 + ]^ 




-13.494 9 










^(4) , ^(4) , J4) (4) (4) 




18.571 4 





052 


663 


= 


ln[(Za)-2] (8^Qg + f Q,) + 


^ In a Qg 


24.277 





000 


502 




i(^gi4 + Ql5 + Ql6) 




-2.323 3 


-0 


022 


511 


(7 = 


^2 ^(2) ^^4 ^(4) ^^5 ^(5) ^^2 


J22_ (2,1) 
rriN 




59 


967 


43(10) 



VII. SUMMARY 

We have obtained relativistic, QED and finite nuclear mass corrections to the magnetic shield- 
ing constant in '^He with the uncertainty of 0.1 ppb, which is caused by neglected QED corrections. 



While our nonrelativistic result is in perfect agreement with the previous one by Drake in [11211 . the 
relativistic correction 0.052 663 is in the moderate agreement with result by Vaara and Pyykko 
in Ref. ||5|] 0.04, which is obtained as the difference between DF LR and HF values in their Ta- 
ble I. We think therefore, that our calculation needs more accurate confirmation, since we have 
not been able to test individual relativistic corrections. However, under the assumption that the 
present calculation is correct, the shielding factor for '^He is now known with the highest accuracy 
of any atom, which supports its use as a NMR standard. Moreover, the theoretical accuracy can 
be further improved by the complete calculation of the QED effects, which for He is certainly 
possible. Regarding calculations for other light atoms and molecules, we are not convinced that 
the commonly used Gaussian functions can be applied for the accurate evaluation of second order 
matrix elements, especially for Qy, Q12 and Q13. It is likely that the use of linear terms, which 
improve the cusp condition of the nonrelativistic wave function, will be necessary in order to ob- 
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tain numerical result with predictable uncertainty. If this can be achieved, it will open a window 
for high accuracy determination of nuclear magnetic moments. Since the NMR frequencies can 
be measured very accurately, as that in ^He, with respect to the proton in tetramethylsilane (TMS) 
ll^, the magnetic moment of helion can be related to the accurately measured proton magnetic 
moment il4Q. 
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